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REDUCTION  THEORY  AND  EXPERIMENTAL 
INJECTED  POLYMER  SOLUTION  UNDER 


INVESTIGATION  OF  IN  SLIT 
HIGH  REYNOLD  NUMBERS 


Wang  Xiliang  Xia  Changsheng 


The  addition  of  a  small  amount  of  drag-reducing  polymer 
in  near-wall  flow  may  greatly  reduce  the  sliding  friction.  In 
engineering  we  are  interested  in  the  theory  and  experimentation 
of  high  Reynold  number  slit  injection  of  drag-reducing  polymer. 
Our  results  are  presented  below. 

1.  Differential  relations  of  wall  shear  stress  of  slit- 
injected  drag-reducing  polymer  solution  in  flow  around  a  stream¬ 
lined  body  of  revolution. 


This  method  is  based  on  the  velocity  distribution  of  the 
joint  boundary  layer  in  references  (4)  and  (5).  The  equation 

is 


M+  «  «-  /liny*  +  B  +  bay*  +  AS 

V* 


(1) 


where  u  is  the  velocity  component  along  x  direction  In  the 

boundary  layer;  y+ -  pmy/v;  p,  -  VtTTp  ls  the  shear 

velocity,  t  is  the  wall  shear  stress,  p  is  the  density;  A, 

™  v  dp 

B  are  constants,  with  A  =  2.5  and  B  =  5.5;  *“■ - ,* 

dx 

v  is  the  kinematic  viscosity,  p  is  the  pressure; 

AS  —  Fin T  —  r (<■)  is  related  to  the  concentration 

c,  v«cr  is  the  shear  velocity  at  the  onset  of  brag-reduction.:. 


for  zero  or  negative  pressure 
gradient 

for  positive  pressure  gradient. 
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The  continuity  equation  of  the  flow  around  a  body  of 
revolution  is 


Q(r»)  ,  d(rt')  — 

6x  dy  (2) 


When  normal  stress  is  neglected,  the  equation  of  motion  is 


ud±.  +  vdlL„-±.d£_+±QLlll 

dx  dy  P  dx  pr  dy 


(3) 


where  x  and  y  are  respectively  the  tangential  and  normal  direc¬ 
tion  of  the  object  profile  in  the  meridian  plane;  r  =  rQ  + 
ycos<|),  4)  is  the  angle  between  the  x  direction  and  the  axis, 

r  is  the  distance  from  a  point  on  the  object  surface  to  the 
o 

axis.  (Figure  1). 


Figure  1.  Coordinates  of  the  body  of 


We  assume 
with  thickness 


the  boundary  layer  to  be  a 

6<<r  . 
o 


Integrating  equation  (2),  we  get 


revolution . 
thin  boundary  layer 


[v,(A\ny*  +  bay*  +  B  ■+•  AZJ))  “  (.u*  +  A  +  bay*') 

dx  dx 

r.  *  <*«  ,  </cab)i  . 


dx  dx 


Substituting  equations  ( h )  and  (5)  into  (3)  and  noticing 


that 


i<*  V  <*y  / 


then  we  get 


+>  ^  ,  *  [  +  y+  _  J-  y«->  - \-^L.v\y*  (u*  —  ~  A')^ 

U  V*  if  +  b  *17  l  y  2  ^  Jy^J  t  dx  *  *  \  2  '  4y 

+du*\d(&B) _ 1  -x.  J_  Kill  -  r;  4-  ^ 

+  'l{**-y*j?)-4T-  -7dx+,r  dy  Udx  +  ^dy* 


Integrating  the  above  equation  from  0  to  6  =  v#6/v  we  get 

iL*G  +  b„l  *!L  h  _  -L  i!L  v\s*  \(a+  ^ \  b'a'S*' 

dx  dx  r  dx  !■'  2/3 

-  i-  AboS*  -  2^1  +  v%  6*  («r  -  baS *  -  2 

-vJTs*-Tr-  (6) 

dx  h 

where  U  is  the  flow  velocity  at  the  outer  edge  of  the  boundary 
layer;  y  is  the  viscosity; 

♦I  -  V 

o  —  «  .  •*  — 

C  -  j -  0*  [(<r  -  A)'  +  A'  -  ^  (<r  -  j  ^)] 

H  -  -  ^,  +  2 j|  V 

Let  the  stiffness  approach  infinity  in  equation  (6)  and  let 
V  •-  V/Vm,  x  —  x/L,  r  —  r/L  ,  where  is  the  flow  velocity  at 

infinity.  L  is  the  length  of  the  body.  Finally  we  obtain  the 
differential  relation  for  the  shearing  stress  of  the  body  of 
revolution  with  slit  injection. 
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(G  -  3ba}»  a'  +  £  o^B*  -  C)  - 

+  j  aS*  [(/<  +  a  -  ~  b*  a1  8*‘  -  -|-  Abas'*  -  2/i'] 

-  aS*  (a  -  2A  —  baS*)  (AB)'  -  RL  V 

(7) 

in  which  the  *  and  the  "  represent  the  first  and  second  order 
derivative  wrt  x; 

Rt  -  UmL/ V  ; 


U»)'-£ont>  +  r(£-£) 


2.  Discussion  on  the  differential  relation  of  the 
shearing  stress. 

1)  when  B=0,  equation  (7)  is  the  shearing  stress  differ¬ 
ential  relation  for  the  flow  around  the  body  of  revolution 
without  the  additive. 

2)  From  equation  (8),  we  have  ( ab —  T  a'i/o,  in  plane 

uniform  solution  flow.  Equation  (7)  becomes 

S*Ua-  Ay  +  A'  4-  r  (a  -  2A)]  a  -  RL  V  or  s*  l(a  -  A )* 

+  ^ -4- r  (<7  -  2^)] -g- -  1 


Equation  (9)  has  the  solution 

R,  -  (o-„rS  [<*(£+1)  +  <r  +  2/1>a(I+i)  ~2A(-r~  A’>ai]  (  1 

when  T/A  is  an  integer.  f 

•  r  •  p 

Here  a,  -  At?,  at_  -  /<  (<m  el  -  ) 

When  r/A  =  0,  this  is  then  the  solution  for  the  flow  around 
the  body  without  the  additional  plane.  q1  plane  d i1  n f 
coefficient  (2)  is  Cf  _  J_['L  iE*  _  JL  [JL  j*+  „♦  (<r  _  .♦)  4 ) 
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For  our  model,  we  have 


r  -  2  A)  (H) 

Rl  a 

Figure  2  shows  the  relational  curves  of  C ^  to  the  Reynold 
number  Re  for  poly ( ethylene  oxide)  solution  of  various  concen¬ 
trations  T/A  =  1,  2,  5  and  for  pure  water  T/ A  =  0.  Here 
=  10  m/s,  V*  =  0.023  m/s. 

3)  For  the  case  of  pure  water,  the  separation  condition 
for  the  flow  around  a  body  is  G-3baH=0;  in  the  presence  of  drag 
reducing  polymer,  the  separation  point  condition  is 

G  —  Iba  H  +  rs *  (<r  —  2A  —  bA  S +)  -  0 

The  extra  third  term  is  usually  greater  than  zero.  When  separa 
tion  exists,  we  ^uy  interpret  this  as  moving  the  separation 


Figure  2.  Curve  of  plane  viscosity  for  uniform  solution 

point  further  to  the  rear  (o,6+,  G  and  H  also  change  with 
additive).  In  practice,  the  curve  for  wall  shear  stress 
coefficient  rmfpVl  —  1/tr1  is  calculated  from  equation  (7). 
In  the  case  of  pure  water,  the  tail  of  the  body  of  revolution 


5 


curve  decrease  since 


2 

makes  o'  increase  rnnidly  and  "i  /p’J 

—  w  2 
V*/V  and  r’/r  are  both  larr.e  negative  numbers.  When  Tv,/pU-*0, 

separation  occurs.  With  additive,  because  of  the  extra  term 

ra+(<r  —  2A  —  bad*)  s  the  coefficient  of  o'  increases  so  that 

the  fall-off  phenomenon  of  the  curve  due  to  V*  /V  and  r'/r  is 

moderated,  delaying  the  occurence  of  separation. 

4)  It  is  very  important  in  solving  equation  (7)  to  estab¬ 
lish  the  wall  concentration  relation  along  the  longitudinal 
direction  and  to  derive  from  which  the  expression  for  (AD)’. 
When  we  assume  that  the  diffusion  of  the  slit  injected  polymer 
at  the  boundary  layer  is  similar  to  the  diffusion  law  of  the 
line  source  at  the  turbulent  boundary  layer,  the  final  stage 
wall  concentration 


cm  GO  -  f/0.55  su  (12) 

r  7 1 

where  g  is  the  unit  length  injection  rate.  When  f  is 

taken  to  be  directly  proportional  to  Cw,  equation  (12)  may  be 
re-written  as 


Tag*  =  constant 


(13) 


From  (13)  with  the  consideration  of  the  expressions  for  a  and 
AB,  we  can  derive 


where 


(AB)'  -  [K,  (<r  +  r  4-  AB  -  3 bad")  -  (r  +  AB)) 

c 

+  ir  -  k,  cr  -  ib *  s+)j  £  -  k,  bas+  —■ 


K,  ~ 


A  B 


A  B  —  A  —  baB* 


(14) 


3.  Experiment  on  drag-reduction  and  wall  concentration 
measurement  for  slit  injection  on  a  plane  and  on  a  body  of 
revolution . 
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1)  The  measuring  procedure 

The  model  used  In  this  experiment  consists  of  a  plane 
and  a  body  of  revolution.  The  length  of  the  plane  is  3.  h  m. 
and  its  thickness  is  25  rrun ,  with  both  ends  trimmed  sharp.  There 
is  a  slit  of  width  0.8  mm  at  an  angle  of  about  30°  with  the 
surface  from  the  front  end  of  the  plane.  During  the  experiment 
the  part  of  the  plane  under  the  water  line  is  0.338  m.  -’ix 
sampling  slits  of  width  0.25  mm  and  effective  length  30  -  35 
mm  are  placed  along  the  longitudinal  direction.  Their  posi¬ 
tions  are  tabulated  in  tables  1  and  2.  The  body  of  revolution 
has  a  length  of  5.463  m.  There  is  an  annular  injection  slit 
of  width  about  0.8  mm  and  at  an  angle  of  about  7°  with  the 
surface  at  a  distance  of  0.439  m  from  the  front  end  of  the 
body.  Seven  sampling  slits  similar  to  those  for  the  plane  are 
placed  along  the  longitudinal  direction.  We  used  the  vacuum 

a 

suction  sampling  method  with  a  sampling  rate  of  about  0.3  cm.  / 
s.crn.  For  concentration  determination  we  used  the  turpidity 
method:  1  cc  50%  sulphuric  acid  and  1  cc  1"  phosphorous  moly¬ 

bdenum  acid  are  added  to  10  cc  of  the  solution  and  the  mixture 
is  shaken  until  the  mixing  is  uniform.  10  minutes  is  allowed 
for  the  solution  to  settle  and  then  its  optical  density  is 
measured  with  a  Model  72  photometer.  The  drag-reducing  agent 
is  poly (ethylene  oxide)  with  a  molecular  weight  of  about  3 
million.  The  injection  concentration  is  approximately  500  ppm 
and  the  injection  rates  for  the  plane  and  the  body  of  revol¬ 
ution  respectively  are  1  kg/s  and  0.74  kg/s. 

2.  Experimental  results 

Figures  3  and  4  represent  respectively  the  experimental 
results  of  the  slit  injection  for  the  plane  and  for  the  body 
of  revolution.  The  experimental  Reynold  numbers  are  respectively 
3  x  10^  -  2.7  x  10 ^  and  5.5  x  10^  -  5.4  x  10^.  The  largest 
reduction  of  drag  coefficient  is  approximately  45£  and  39% • 

For  constant  injection  rate  and  injection  concentration,  the 
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Figure  3*  Drag  experiment  Figure  .  Drag  experiment 

curves  of  the  plane  for  pure  curve  of  a  body  c-f  revolution 

water  and  for  slit  injection  for  slit  injection  of  poly 

of  500  ppm  poly (ethylene  (ethylene  oxide)  solution  and 

oxide).  the  comparison  between  theo¬ 

retically  calculated  percentag 
drag  reduction  and  the  experi¬ 
mental  results. 


drag  reduction  decreases  slightly  with  increasing  velocity. 

This  decrease  is  more  pronounced  for  the  plane  than  for  the 
body  of  revolution  because  at  high  velocity  there  is  less  water 
immergence  for  the  plane  due  to  the  free  surface  effect. 


The  measured  results  of  the  wall  concentration  for  the 

plane  and  for  the  body  of  revolution  with  slit  injection  are 

tabulated  in  Table  1  and  Table  2  in  which  V  is  the  flow  speed 

(m/s),  Q.  is  the  injection  rate  (kg/s),  C  is  the  injection 
J  J 

concentration  (ppm),  Cw  is  the  wall  concentration  (ppm),  and 

x  is  the  distance  (m)  from  the  front  end.  *  indicates  a  higher 

concentration  because  the  sampling  bottle  was  not  thoroughly 

cleaned  after  sampling  for  the  Cj  =  1000  ppm  case. 

The  ^results  for  the  wall  concentration  shows  that  the 
variation  of  wall  concentration  is  larger  near  the  slit  outlet. 
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The  wall  concentration  within  1.5  n;  fro::.  the  .oil1  is  bi,:-:er 
than  the  saturated  concentration  of  25  ppm  for  drar  redact  io: 
There  seers  to  be  a  re 'ion  with  relatively  so  all  variation 
in  the  concentration  for  a  fairly  lenr  distance  toward  the 
rear,  especially  for  the  case  of  the  body  of  revolution.  It 
is  difficult  to  reduce  the  wall  concentration  to  an  express¬ 
ion  used  in  reference  [3],  especially  for  the  rc-'ion  where 
C.  =  l.  For  the  case  as  described  in  our  paper,  the  wall  ccn- 
centration  must  be  affected  by  the  injection  ancle,  the 
pressure  and  the  surface  curvature  of  the  object. 

Table  1  Wall  Concentration  for  the  plane  with  injection 

—  conditions  , . .  ■■  — — ■  -- —  ■  —  — — — ; — —  - ;■■■  1  - 


V- 3.540 

V  =  5.106 

V  -  6.561 

V  =  7.620 

V  -  3.480 

£>,  -  1.0 

£>,-1.0 

£>,  =  1.0 

0,  =  1.0 

£>,-1.0 

C,  -  536 

C,  =  536 

C,  =  536 

C,  =  536 

C,  -  1000 

Table  2  Wall  concentration  for  the  body  of  revolution 

with  injection. 


conditions 
of  injectio 

vc.  \ 


n  V  =  4 . 060 

V  «=  5.SS9 

V  ~  6.962 

V  —  4.230 

£>,  =  0.735 

£>,  =  0.735 

£>,  »  0.735 

£>,  -  0.735 

C,  =  536 

C,  =  588 

C,  -  588 

C,  -  1000 

104 

200 

233 

793 

116 

112 

100 

530 

114 

104 

110 

246 

•  12.6 

5.9 

7.7 

118 

1.1 

5.5 

6.7 

14.7 

6.1 

5.4 

6.2 

8.5 

6.3 

6.0 

6.2 

8.5 

t'  -  4.055 
£>,  •»  0.89 
C,  -  588 


4)  Comparison  between  theoretical  calculation  and 
experl::. ent: al  result 

We  carried  out  theoretical  calculation  for  the  body 
of  revolution  of  length  5.463  m  at  velocities  =  4,  €  and  8  m/s. 
The  method  in  reference  [8]  was  used  for  potential  flow  cal¬ 
culation.  Runre-Kutta  method  was  used  to  solve  quation  (7). 

When  we  took  C  =m/0.555U  and  used  enuation  (1J;)  for  ( A B )  ’  as 
w  ' 

derived  from  it,  the  value  for  the  drag  reduction  was  somewhat 

lower  than  the  actual  value.  In  this  paper  we  used  in  our 

calculation  the  average  value  of  the  wall  concentration  as 

determined  by  the  experiment.  The  initial  value  was  determined 

from  °  —  wina'’’  +  bae*  +  B  +  ab  where  the  boundary  layer 

thickness  5  «,  q  059S//(iogRes  —  3.17)  ,  s  is  the  arc  length,  and 

( 9 ) 

the  corresponding  Reynold  number  is  Res  .  In  the  calcula¬ 
tion,  Afl  —  C.log  -^2-  .  VJhen  the  velocity  was  4m/s,  the 

curves  for  the  shear  stress  coefficient  r ^/pU1  —  l/<^  are 
shown  in  Figure  5  for  pure  water  and  for  poly (ethylene  oxide) 
with  injection  concentration  500  ppm. 


2 

Figure  5-  Theoretically  calculated  curves  of  xw/pU  for 
water  and  for  solution.  ^ 

Integrating  along  the  surface  for  wall  shear  stress  for 
water  and  for  solution  as  well  as  for  their  difference,  we 
obtain  respectively  the  value  for  drag  coefficient  and  drag 
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reduction.  We  also  plotted  the  theoretically  calculated  per¬ 
centage  drag  reduction  In  Figure  from  v;hlch  one  can  see  that 
the  theoretical  calculation  basically  agrees  with  experimental 
result.  This  method  may  be  used  to  predict  the  calculated 
drag  reduction  for  a  body  of  revolution  with  injection. 
Strictly  speaking,  when  applying  Meyer's  model  to  external 
flow,  and  in  particular  to  external  flow  with  uniform  concen¬ 
tration,  the  assumed  form  AB  —  r(C„)ln  -£*-  requires  proof 

fmu 


by  direct  measurement. 


The  terms  in  equation  (7)  of  this  paper  have  clear  physi¬ 
cal  meaning  and  its  calculation  is  simple.  Therefore  it  fur¬ 
nished  us  a  method  to  calculate  drag  reduction  that  basically 
agrees  with  experimental  result. 


We  are  obliged  to  Comrades  Zhou  Dexiang  and  'Wang  Gueiqin 
for  supplying  data  from  concentration  measurement. 
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Notation:  rc ,  R,f)t  A*nRr  /ifcf^  denotes  respectively  the 


outer  radius  of  the  centrifugal  nozzle  outlet,  the  vortex  flew  (- 

arm  through  which  the  fuel  enters  the  vortex  flow  chamber,  the  * 

inlet  cross-sectional  area,  and  the  geometrical  characteris- 

p 

tic  parameter;  f. '  ,  w,  a,  ,  <->=1- ( rv/rg )  denote  respectively  j 

the  discharge  rate  of  the  simple  centrifugal  nozzle,  the  dis-  j 

charge  coefficient,  the  outlet  spray  angle,  the  radius  of  air  ; 

vortex,  and  the  coefficient  of  effective  outlet  cross-section; 

o„  <Pi  “  i  ■"  denote  respectively  the  discharge  rate  I 

of  a  centrifugal  nozzle  with  annular  cross-sectional  outlet,  3 

the  discharge  rate  coefficient,  the  outlet  spray  angle,  the 

radius  of  the  central  cylinder,  and  the  coefficient  of  annular  'j 

outlet  cross-section  area;  w,  w  ,  wm,  w  denote  respectively 

a  T*  o  1 

the  total  velocity  at  an  arbitrary  point  r  in  the  outlet  cross- 

sectional  area,  the  axial  velocity,  the  tangential  velocity  and 

the  total  velocity  when  the  pressure  energy  is  completely  con-  , 

verted  into  kinetic  energy;  p  ,  p  ,  p’  denote  respectively  the 

oil  pressure  in  front  of  the  nozzle,  the  static  pressure  at 

the  point  where  the  centrifugal  nozzle  with  annular  outlet 

r  =  r  ,  and  the  static  pressure  at  the  point  where  the  simple 
s 

centrifugal  nozzle  outlet  r  =  r  (all  are  residual  pressures 

s 

relative  to  the  ambient  medium);  \p  is  the  corrected  discharge 

coefficient  of  the  nozzle  inlet  passage  and  the  coefficient 

of  flow  column  shape  change  as  the  fuel  enters  into  the  vortex 

chamber,  y  is  the  fuel  weight,  and  g  is  the  gravitational 

acceleration;  the  subscripts  1,  2,  3  denote  respectively  the 

results  corresponding  to  the  3  different  assumptions  on  p  , 

s 

s  denotes  the  value  at  the  point  where  r  =  r  or  denotes  the 

s 

value  of  the  annular  outlet  cross-sectional  area,  the  sub- 
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the 


I 


script  v  denotes  the  value  at  the  air  vortex 
superscript  denotes  the  vniue  corresponding  to  the  case  when  1 

the  central  cylinder  is  absent  (simple  centrifugal  nozzle).  | 


1.  Centrifugal  nozzels  are  widely  used  in  power,  aero¬ 
nautic,  chemical  and  metallurgical  industries.  The  operational 
characteristics  of  the  simple  centrifugal  nozzle  has  Ions  been  , 

analysed  by  G.  ’!.  Abramovich  [1  ]  (we  shall  refer  to  it  as  A's 
theory).  A  large  amount  of  experimental  and  theoretical  mat¬ 
erials  from  analysis  have  been  accumulated.  Ordinarily  for  j 

furnace  nozzles  the  experimental  results  do  not  deviate  from 
A's  theory  significantly.  This  shows  that  the  theory  basically  j 

reflects  the  operational  characteristics  of  the  centrifugal  ( 

nozzle.  Therefore  it  remains  to  be  the  foundation  of  indus-  I 

trial  calculations,  and  forms  the  point  of  departure  for  the 

[2] 

calculation  of  more  complicated  nozzles  with  structures.  t 

I 

To  increase  the  range  of  centrifugal  nozzle  flow  rate 
regulation,  various  structures  involving  such  forms  as  fuel 
return,  nozzle  outlet  cross-sectional  area  regulation,  or  co-  j 

axial  placement  of  principal  and  auxilliary  nozzles,  etc.  may 
be  used  [2,3]  .  In  several  of  these  cases,  we  encounter  the 
centrifugal  nozzle  with  annular  outlet  cross-section .  Research 
on  this  kind  of  nozzles  at  present  is  still  rather  scanty. 

In  this  paper  we  shall  analyse  a  nozzle  which  relies  on  a 
change  in  the  cross-sectional  area  of  its  annular  outlet  to 
regulate  the  fuel  rate  (Figure  1).  Its  only  difference  from 
the  simple  centrifugal  nozzle  is  in  the  installation  of  a 
conical-shaped  pin  rod  at  the  center.  The  fue*l  rate  is  regu¬ 
lated  by  moving  the  pin  rod  forward  or  backward  so  as  to  change 
the  size  of  the  annular  cross-sectional  area  of  the  outlet 
nozzle,  while  keeping  the  fuel  pressure  constant.  The  typical 
result  of  experiments  performed  on  the  nozzle  test  platform 
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Figure  1.  Centrifugal  nozzle 
with  adjustable  outlet  cross- 
sectional  area. 

1.  Orifice,  2.  Vortex  flow 
disc,  3.  Centra]  pin  rod 
(can  be  moved  forward  or 
backward),  fuel-separation 
ring. 


Figure  2.  Regulation  char¬ 
acteristics  of  centrifugal 
nozzle  with  adjustable  outlet 
cross-sectional  area. 

Rectangular  tangential  slots 
of  dimension  1.87  x  2.595 
A  *  1.575 


is  shown  in  Figure  2.  From  Figure  2  we  can  see  that  when  the 
pin  rod  is  being  moved  forward  and  the  blockage  of  the  outlet 
is  being  gradually  increased,  there  is  an  initial  stage  during 
which  the  nozzle  flow  rate  and  the  spray  angle  are  practically 
unaffected  by  the  blockage.  They  only  show  significant  change 
when  the  degree  of  blockage  is  greater  than  a  certain  value: 
a  rapid  decrease  with  increasing  blockage.  Similar  experi¬ 
mental  result  is  obtained  in  reference  [3]. 

2.  For  conciseness,  we  shall  base  our  analysis  in  this 
paper  on  the  theory  of  an  ideal  fluid  in  a  simple  centrifugal 
nozzle.  However,  in  order  to  obtain  better  accuracy  in  our 
calculation  with  the  theoretical  formulae,  we  shall  first  make 
some  modifications  to  the  original  A's  theory  on  simple  cen¬ 
trifugal  nozzles.  ' 

In  a  simple  centrifugal  nozzle,  owing  to  the  centrifuge 
effect  produced  by  the  tangential  entry  of  the  fuel  oil,  the 
fuel  does  not  flow  through  every  part  of  the  nozzle  outlet 


| 


ci’oss-uect  ionai  area .  Instead ,  an  air  vortex  is  formed  near 
the  axis  of  Uic  nozzle,  with  the  coefficient  of  effective 
c  ross -s  oe  t i c n a  1  area  <p  “  l-Oa/O’^l  .  The  magnitude  of  <p 
is  Jo t or;:'. irou  by  A,  the  geometrical  characteristic  parameter 
of  the  nozzle  at  infinite  stiffness: 


A 


x  R  r{ 

4>fi 


(1) 


The  geometrical  characteristic  parameter  A  of  the  nozzle, 
the  coefficient  of  effective  outlet  cross-section  area  <J> ,  the 
nozzle  discharge  coefficient  and  the  spray  angle  a  are  related 
as  follows: 


A 


b 


a 


2  sin 


(2) 

(3) 

(4) 


We  differ  from  the  original  A’s  theory  on  only  two  counts:  one 
is  that  we  introduced  in  equation  (1)  the  correction  coeffici¬ 
ent  to  correct  the  effect  introduced  by  the  fact  that  the  dis¬ 
charge  coefficient  in  the  tangential  entrant  passage  is  not 
equal  to  1  and  the  fact  that  the  fuel  shape  deforms  as  the 
fuel  enters  the  vortex  flow  chamber,  two  is  that  we 

have  taken  into  consideration  the  effect  on  the  spray  angle 
when  the  static  pressure  of  the  fuel  is  transformed  into  axial 
velocity  at  the  nozzle  outlet.  Let  sin(a/2)  =  (average  tangen¬ 
tial  velocity  at  the  nozzle  outlet )/(the  total  velocity  when 
the  pressure  before  the  nozzle  is  completely  transformed  into 
kinetic  energy )  . 


Table  1  compares  the  experimental  results  with  the  results 
calculated  from  equations  (1)  -  (4).  The  average  value  and 
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upper  and  lower  error  Units  for  16  nozzles  are 

•£**■  -  1.004  -  1.015  !?;!! 

while  the  ratio  of  the  experimental  results  with  that  from  the 
original  A's  theory  are 

-±s*-  -  0.88  «;«,  - 0.87  i?;* 

/***•••  ««*■» 

Hence  the  modified  or  corrected  formulae  (1)  -  ( ^ )  agrees 
better  with  experiment  than  the  original  A’s  theory. 

The  nozzles  listed  in  Table  1  all  came  from  actual  produc¬ 
tion  and  therefore  may  be  regarded  as  representative.  For 
the  2  nozzles  used  in  our  experiment,  when  ip  is  taken  to  be 
0.8,  the  calculated  results  from  (1)  -  (A)  also  agreed  well 
with  experiment.  When  the  central  pin  rod  is  ineffective 
(equivalent  to  a  simple  centrifugal  nozzle),  for  the  nozzle  with 

3  tangential  circular  holes  and  A  -  1.123,  the  discharge 
coefficient  has  a  computational  value  of  0.^11,  the  experimen¬ 
tal  value  being  0.371,  computed  value  for  the  spray  angle  is 
69.5°,  the  experimental  value  being  61°;  For  the  nozzle  with 

4  rectangular  tangential  slot  and  A  =  1.575,  the  computed  value 
for  discharge  coefficient  is  0.337,  the  experimental  value 
being  0.335,  the  computed  value  for  the  spray  angle  is  78.7°, 
the  experimental  value  being  72°. 

In  Table  1,  we  take  =  t.8  to  compute  A.  Strictly 
speaking,  ip  should  be  dependent  upon  the  structure  of  the 
nozzle  and  its  state  of  mechanical  processing.  For  simplicity, 
we  shall  take  ip  =  0.8  in  all  the  subsequent  analysis. 

For  computation  of  equations  (2)  and  (3)  it  is  possible 
to  use  the  curve  in  reference  [1],  Equation  (2)  -  (4)  may  also 


-MU 


Table  1. 


Key:  (1)  Nozzle  No.;  (2)  Average;  (3)  Geometric  Character¬ 
istics;  (4)  Discharge  rate;  (5)  Ratio  of  experimental 
to  theoretical;  (6)  Spray  angle. 


be  approximated  with  the  following  relations: 


0.428  A-*  M 

(5) 

0.630 

(6) 

65.2°  A 0  M7 

(7) 

We  obtained  these  approximate  relations  by  using  the  least 
square  method  on  a  numerical  calculator.  In  the  region 
A=0.7  -  4.0,  the  results  do  not  differ  from  those  obtained  by 
using  equations  (2)  -  (4)  by  more  than  2%. 


3.  Let  us  now  analyse  the  operational  characteristics 

of  the  centrifugal  nozzle  with  an  annular  outlet.  Analogous 

to  the  treatment  in  reference  [1],  we  consider  the  case  of 

an  ideal  fluid.  Let  the  radius  of  the  nozzle  outlet  be  r  , 

c 

and  the  radius  of  the  concentric  cylinder  be  rg  so  that  the 

actual  nozzle  outlet  cross-sectional  area  is  the  annular  region 

bounded  by  the  concentric  circles  with  radii  r  and  r  (Figure 

c  s 

3).  There  are  two  possible  cases:  (1)  r  is  less  than  the 

s 

radius  of  the  air  vortex  as  calculated  from  the  A  value;  (2) 
rg  is  larger  than  this  vortex  radius.  We  discuss  the  cases 
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separately  below: 

1)  Case  of  ry  (Figure  3a).  Here  the  cylinder  only 

occupies  a  part  or  the  whole  of  the  vortex  cross-sectional 
area,  and  will  not  affect  the  flow  of  the  oil.  Therefore  the 
discharge  coefficient  and  the  spray  angle  of  the  nozzle  are 
exactly  the  same  as  in  the  absence  of  the  cylinder.  VJe  can  use 
the  formulae  (2)  -  (4)  for  the  simple  centrifugal  nozzle  in 
our  computation.  This  prediction  has  been  verified  by  our 
experiment:  from  Figure  2,  when  (r,/r,y<0.38  ,  the  dis¬ 

charge  coefficient  and  the  spray  angle  of  the  nozzle  are 
independent  of  (r  /r  )2  when  the  computed  value  of  (r  /r  )2 
is  0.384,  which  demonstrates  that  for  rg  <  r  ,  the  magnitude 
of  rg  does  not  affect  the  discharge  coefficient  and  the  spray 
angle  of  the  nozzle.  The  experimental  result  of  reference  [3] 
agrees  with  this  conclusion. 

2)  Case  of  rg>rv  (Figure  3b).  Here  the  insertion  of  the 
central  cylinder  not  only  will  make  the  air  vortex  disappear, 


Figure  3.  Outlet  Cross-section  of  a  centrifugal  nozzle 
with  an  annular  outlet 
(a)  rs<rv  (b)  rg  >  ry 


but  will  also  directly  affect  the  area  of  the  oil  flow.  Let 
G  and  respectively  represent  the  flew  rate  and  the  discharge 


coefficient-  of  the  c^ntr 1  t’ugai  no-^lo  with  the  annular  outlet 

p  be  the  residual  pressure  of  the  fuel  ell  in  the  nozzle 
o 

relative  to  the  ambient  .'tedium,  y  be  the  weight:  of  the  fuel 
oil,  arid  £  be  the  gravitational  acceleration,  then 


C. 


h‘*rW  2 


nrp» 


(8) 


It  is  easy  to  prove  with  the  same  method  as  in  reference  [1] 
that:  the  axial  velocity  in  the  cross-sectional  area  of 

the  outlet  is  still  uniformly  distributed,  independent  of  r. 

G, 

Thus,  wt  — - ; — ,  where 


y  uioi-i  j-  uu  t' 


r*  r.<p, 

We  notice  that  the  total  velocity  at  the  inside  boundary  r 

s 

of  the  nozzle's  annualr  outlet  (assuming  Ps  to  be  the  resi¬ 
dual  static  pressure  with  respect  to  the  ambience  at  that 
point)  is 


in  which  the  axial  velocity 
tial  velocity  - - OlA - G> 


and  tangen- 


Wj,  — 


rfi4>r,  V*r\J  1  — 


Hence 


I -  Jl 

G,  —  r*  T\  o  — t  ==^=j= 

V 


(9) 


Comparing  (8)  and  (9),  we  know  that  the  discharge  coefficient 
of  the  centrifugal  nozzle  with  annular  outlet  when  r  >r  is 

S  V 


»-'F-f.Nh  +  :t=7. 


For  a  given  nozzle  set  of  and  operating  conditions,  p  ,  A,  <)> 
are  all  known  quantities.  Therefore  the  key  to  computing  ug 
is  the  determination  of  the  residual  pressure  p  at  the  surface 

o 

r  .  Generally  speaking,  p  may  be  a  complicated  function  of 
s  s 

<J>  (or  r  ).  Since  at  present  there  is  a  lack  in  direct  experi- 
s  s 

mental  data  about  p  ,  we  adopt  the  method  of  making  all  kinds 

s 

of  assumptions  on  p  ,  deriving  the  corresponding  formula  for 

s 


T 
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U  and  then  comma ring  the  different  commuted  results  with  the 
experiment  to  determine  which  assumption  about  p  is  more  in 
agreement  with  reality. 


We  have  made  the  following  three  assumptions  about  pg: 

(1)  We  assume  that  the  static  pressure  at  r=rc^  stays  the  same 
before  and  after  int reducing  the  central  cylinder.  Thus  the 
static  pressure  po^  at  tho  surface  of  the  cylinder  is  equal  to 
p’,  the  static  pressure  at  r=r  in  a  simple  centrifugal  nozzle. 
Analogous  to  reference  [1] 

_  i  -  „?.)  -  i (£*Y(i  _  X)  _  <„ (_> - _J_) 

2g  2g  'Y'l’fi'  'rti  rj/  •  \  1  —  V  1  —  <p*' 

whence 


fa 


Ji-to-li-  -  _1_V 

Y  pt  V  v  1  —  V  1  —  9V 

Substituting  into  equation  (10),  we  get 

-  t=t) 

JT73Z 

V  V?  1-9. 


(11) 


(2)  We  assume  that  the  axial  velocity  at  the  nozzle  outlet 

stays  constant  before  and  after  the  introduction  of  the  cylinder, 

i.e.  w  /w’  =  1.  Since  .  .  ,  . 

a  a  JZ*-  -  {<L\  (V\  HdSL 


w,  \Gv  \<p,/  n<p. 


then  we  have 


/‘ai  *“  /‘<P«/<P 


(12) 


Substituting  equation  (12)  into  (10),  we  see  that  this  is 
equivalent  to  assuming 
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0  on  the  surface  of  the 


(3)  We  assume  that  the  pressure  p  0  = 

central  cylinder  at  all  times,  i.e.  even  r  >r  ,  the  static 

s  v 

pressure  on  the  cylindrical  surface  is  always  equal  to  the 
ambient  pressure.  From  equation  (10)  we  knew 


(13) 


All  of  the  equations  (11)  -  (13)  derived  from  the  three  differ¬ 
ent  assumptions  give  results  in  qualitative  agreement  with  the 

experiment,  namely  that  when  r  >r  ,  as  r  is  increased  ( d>  is 

s  V  s  s 

decreased),  both  the  discharge  coefficient  and  the  spray  angle 

of  the  centrifugal  nozzle  with  annular  outlet  decrease  in  value 

(see  below).  To  determine  which  assumption  about  p  is  in 

s 

quantitative  agreement  with  the  experiment,  we  compare  the  com¬ 
puted  results  of  equation  (11)  -  (13)  with  experiment.  Figure 
4  and  Figure  5  present  the  ratio  of  the  computed  to  the  experi¬ 
mental  results  for  the  tangential  circular  hole  type  nozzle  and 
the  tangential  rectangular  slot  type  nozzle  respectively. 

From  Figure  4  we  can  see  that  for  this  tangential  hole  type 

nozzle,  when  <f>  is  large  (small  blockage),  the  experimental 
s 

points  follow  more  closely  to  the  curve  calculated  from  equa¬ 
tion  (12)  and  when  <f>  decreases  (blockage  increases),  the  experi¬ 
mental  points  gradually  shift  over  to  the  curve  of  equation  (11). 
From  Figure  5  we  can  see  that  for  the  tangential  slot  type 
nozzle  the  experimental  points  follow  more  closely  the  calculated 
curve  of  equation  (11). 

The  author  of  reference  [3]  changed  <J>  by  installing  con- 

s 

centric  cylinders  of  different  diameters  in  the  center  of  the 
tangential  hole  type  nozzle  when  he  investigated  the  problem 
of  the  Interaction  between  the  principal  and  auxiliary  nozzles 
in  a  gas-fired  turbine.  His  result  showed  that  the  experimental 


points  follow  the  curve  of  equni  icr.  ( ]  2 ) { Fi gure  6 )  . 
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computed  to  experimental  results 


Figure  ^ .  Discharge  coeffi¬ 
cient  of  centrifugal  nozzle 
with  adjustable  outlet  area 
of  the  tangential  hole  type. 

three  #3.2  tangential  circu¬ 
lar  holes  A  =  1.123 


Figure  5.  Discharge  coeffi¬ 
cient  of  centrifugal  nozzle 
with  adjustable  outlet  area 
of  tangential  slot  type. 

four  rectangular  tangential 
slots  1.8?  x  2.595  A  ---  1.575 


Similar  to  the  above,  for  r  >  r  the  spray  angle  is  calculated 

S  V 


from  the  following  equation: 


*(?H 


where  w^  is  the  average  tangential  velocity  at  the  outlet  of 

the  centrifugal  nozzle  with  annular  outlet,  which  we  take  to 

be  the  tangential  velocity  at  (r  +  r  )/2,  i.e. 

s  c 


.  -  £JL  _ _ Idfh _  lit 

1  +  J  1~<p<  '  r 


The  total  velocity  w  when  all  the  pressure  energy  is  trans¬ 
formed  into  kinetic  energy  is 


■  i^ki  ijj  i 


Substitute  the  different  i:c,  derived  from  the  different 
assumptions  (equation  11  -  13)  into  equation  14,  we  get  the 
corresponding  formulae  for  a  : 


u)  rs  =  p,s  *->:■*» 


o.l 


—  2«n~*  (  ■  ) 

n  +  Vi  -  <p/ 


(15) 


(2)  Constant  axial  velocity  before  and  after  introducing 
the  cylinder  ,  . 

„  •  -i  /  2/  ;<,»  \ 

(16> 


(3)  Ps  "0 


o.l 


2  sin"*  ( - 

Vl  +  V 1  -  <p/ 


(17) 


Figure  6.  Ratio  of  discharge  coefficient  of  centrifugal  nozzle 
of  tangential  orifice  type  with  annular  outlet  to  the  computed 
result  from  equation  12. 

four  circular  inlet  holes 


Hence 


2  sin 


-i  7  7Ali'  _ \ 

'1  +  \/l  —  <pj 


(14) 


Substitute  the  different  \i^  derived  from  the  different 

o 

assumptions  (equation  11  -  13)  into  equation  14,  we  get  the 

corresponding  formulae  for  a  : 

s 


(1)  Ps  =  p’s  P.-P.tt>  -  2  sin-*  (  ■  -  - ) 

N1  +  VI  —  <P. 


(15) 


(2)  Constant  axial  velocity  before  and  after  introducing 
the  cylinder 

2  A  f‘,i  \ 


(3)  Ps  =0 


„  —  2  sin"'  ( - f 

'1+vl  -  <P.' 

a„  -  2*in"‘  ( - 

VJ  +  VI  -  W 


(16) 


(17) 


Figure  6.  Ratio  of  discharge  coefficient  of  centrifugal  nozzle 
of  tangential  orifice  type  with  annular  outlet  to  the  computed 
result  from  equation  12. 

four  circular  inlet  holes 
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For  Uk 


ntial  hole  typo  and  the  tangential  slot  type 


no:::: leu  that  we  used  in  our  experiments,  the  ratio  of’  the 
experimental  result  to  that  computed  from  equation  (15)  -  (17) 


are  shown 


ely  in  Figures  7  ant 


8.  It  can  bi 


;een  that 


the  different  formulae  all  give  results  in  qualitative  agreement 

with  experiment ,  na::.ely  that  when  decreases,  the  spray  angle 

a  also  decreases.  Fxeeot  for  small  values  of  4  when  the 
s  ■  s 

error  is  relatively  large,  the  experimental  points  for  the  tan¬ 
gential  hole  type  nozzle  are  closer  to  the  computed  result  from 
equation  (1.6)  while  those  for  the  tangential  slot  type  nozzle 
are  closer  to  the  computed  result  from  equation  (15)  (Figure  8). 
The  principal  reason  why  the  computed  values  deviate  from 
experiment  for  small  8g  is  that  the  parts  of  the  nozzle  are  not 
truly  concentric.  The  nozzles  used  in  our  experiment  are  used 
in  actual  production.  There  is  always  some  degree  of  eccentri¬ 
city  between  the  central  pin  rod  and  the  orifice.  When  the 
blockage  is  large  (or  when  <6  is  small),  the  spray  cone  deviate 
significantly  from  axial  symmetry  so  that  in  some  directions 
the  spray  angle  is  much  larger  than  that  in  some  other  directions. 
This  also  affects  the  discharge  characteristics  of  the  nozzle, 
and  may  well  be  what  caused  the  discrepancy  in  the  results  for 

small  <f>  in  Figure  A  and  6. 
s 


Thus,  for  the  tangential  hold  type  nozzle,  the  discharge 

coefficient  (Figure  A  and  6)  and  the  spray  angle  (Figure  7) 

both  support  the  assumption  that  the  introduction  of  the  central 

cylinder  does  not  affect  the  axial  velocity  at  the  nozzle  outlet; 

however,  for  the  tangential  slot  type  nozzle,  the  discharge 

rate  (Figure  5)  and  the  spray  angle  (Figure  8)  both  support  the 

assumption  that  the  static  pressure  at  r  is  not  affected  by 

s 

the  introduction  of  the  central  cylinder.  Further  investigation 
is  necessary  to  find  out  why  there  is  a  difference  in  the  regula¬ 
tion  characteristics  of  the  two  types  of  nozzles  with  different 
structures . 
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1  is  calculated  Kith  equation  (15);  2  with  equation  (16); 
3  with  equation  (17) 


Figure  7.  The  ratio  of 
calculated  to  experi¬ 
mental  values  of  the  spray 
angle  for  a  centrifugal 
nozzle  of  the  tangential 
hole  typo  with  adjustable 
outlet  cross-sectional 
area . 

Three  3.2.  tangential  holes 


Figure  8.  The  ratio  of 
calculated  to  experimental 
values  of  the  spray  angle 
for  a  centrifugal  nozzle  of 
the  tangential  slot  type  with 
adjustable  outlet  cross- 
sectional  area. 


Four  1.87  x  2.595  rectangular 
tangential  slots. 


4.  In  summary,  the  operational  characteristics  of  the 
centrlfural  nozzle  with  annular  outlet  are  ar  follows: 

1)  When  the  radius  of  the  central  cylinder  is  not  larger 
than  the  value  of  r  computed  from  the  value  of  A  (equation  1), 
the  introduction  of  the  central  cylinder  does  not  affect  either 
the  discharge  rate  or  the  spray  angle  of  the  nozzle.  Computa¬ 
tions  may  be  carried  out  with  the  equations  for  the  simple 
centrifugal  nozzle,  equations  (2  -  4)  (or  c  rial : , ns  (5  &  7)). 
For  the  nozzle  with  co-axially  arranged  principal  and  auxiliary 
nozzles,  the  outer  diameter  of  the  auxllliary  orifice  should  no 
be  larger  than  the  air  vortex  of  the  principal  nozzle  (calcula¬ 
ted  from  equation  2). 

2)  When  r  >r  ,  as  <J>  decreases  (blockage  increases),  the 
discharge  coefficient  and  the  spray  angle  of  the  centrifugal 
nozzle  with  annualr  outlet  both  decrease.  The  decrease  of  the 
discharge  rate  is  necessary  for  regulating  the  fuel  rate,  but 
the  decrease  in  the  spray  angle  is  usually  undesirable.  It  is 
just  because  of  this  that  not  too  many  nozzles  with  only  adjust 
able  nozzle  outlet  cross-sectional  area  are  used  in  practice. 
For  such  nozzles,  the  range  of  the  regulation  ratio  should  not 
be  too  large,  i.e.  <J>  should  not  be  made  too  small,  so  that  too 
small  a  spray  angle  may  cause  an  overly  concentrated  fuel  dis¬ 
tribution  and  poor  atomization.  Within  the  practical  range  of 
regulation  ratio,  we  suggest  that  equation  (12)  and  (36)  be 
used  to  compute  the  discharge  rate  and  the  spray  angle  for  a 
tangential  hole  type  centrifugal  nozzle  with  adjustable  outlet 
cross-sectional  area  and  that  equation  (11)  and  (15)  be  used 
for  the  tangential  slot  type  nozzle. 

If  the  nozzle  inlet  cross-sectional  area  is  also  adjusted 
at  the  same  time  when  the  outlet  cross-sectional  area  is 
adjusted,  then  the  problem  of  the  changing  spray  angle  when 
only  the  outlet  cross-sectional  area  is  adjusted  may  be  over- 
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come.  If  d-.-s  li.n-.-J  j. -rlv ,  the  spray  an.:!--  may  wo  1  1  remain 
unchanged  win-n  *-1/'  discharge  rate  of  t  he  nozzle  Is  In*  inn 
adjusted .  Some  e xr  erl  mental  result  in  this  direction  was  given 
in  reference  [•'*],  inJicatinr  that  this  is  really  an  improvement 
with  practical  value,  h  w-  v-  r,  i--. ’'-r-n.--  [hj  failed  to  give  the 
correct  ec:: ; u\  at!  onnl  so  heme  for-  this  case  of  siniultaneously 
adjusting  both  the  inlet  and  outlet  cross-sectional  areas.  A 
slight  modification  in  the  computational  scheme  in  our  paper 
should  enable  us  to  apply  it  to  this  more  complicated  situa¬ 
tion;  we  need  only  to  consider  that  in  the  regulating  process, 
as  varies,  the  inlet  cross-sectional  area  also  varies,  hence 
the  geometrical  characteristic  parameter  A  calculated  from 
equation  (1)  is  now  a  variable.  The  effect  of  this  computational 
scheme  will  be  left  for  future  investigation. 

Addendum:  The  cool  state  experimental  result  of  the  power 

plant  fuel  nozzle  cited  in  Table  1  was  obtained  jointly  by  the 
Beijing  Central  Research  Center,  Department  of  Hydroelectricity, 
the  Beijing  Thermal  Power  Generating  Station  and  the  Xian  Ther¬ 
mal  Engineering  Research  Center.  We  also  cited  part  of  the 
cool  state  experimental  result  on  nozzles  with  adjustable  out¬ 
let  cross-sectional  area  done  at  the  Capital  Steel  Forging 
Plant.  Comrades  Jin  Yongli,  Zhang  Buzhou,  Zui  Xiansheng,  Zhang 
Kerning,  Huang  Jiaqi,  Yu  Zhemin  and  Wang  Shengli  also  partici¬ 
pated  in  this  experiment. 
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